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Abstract. The main steps involved in realistic shell- model calculations employing 
two-body low-momentum interactions are briefly reviewed. The practical value of this 
approach is exemplified by the results of recent calculations and some remaining open 
questions and directions for future research are discussed. 
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1. Introduction 

A key challenge for nuclear physics is to develop a comprehensive understanding of the 
observed properties of nuclei starting from the forces among nucleons. Nowadays there 
are two main ways to try to achieve this goal, which are substantially complementary 
to each other. The first one is based on the so-called ab initio calculations in which 
nuclear properties, such as binding and excitation energies, are calculated directly from 
first principles of quantum mechanics using an appropriate computational scheme. To 
this category belong the Green's function Montecarlo Method (GFMC) fl], the no- 
core shell- model (NCSM) [2], and the coupled cluster method |3]. Clearly, all ab initio 
calculations need huge amount of computational resources and are therefore currently 
limited to light nuclei. 

The second line of attack consists of using the shell model with realistic two-body 
effective interactions, namely derived microscopically from the free nucleon-nucleon 
(NN) potential. In this kind of calculations only the valence nucleons are treated as 
active particles while the core polarization effects are taken into account perturbatively 
in the derivation of the effective interaction Ves- Of course, this approach allows one to 
perform calculations for medium- and heavy-mass nuclei which are far beyond the reach 
of ab initio calculations. 

In this context, it should be mentioned that in recent years successful shell- model 
calculations, in particular for sd- and p/-nuclei, have been performed [H El [6] by 
employing modified versions of realistic effective interactions. Let us make it clear that 
this paper is concerned only with calculations employing two-body effective interactions 
derived from the free NN potential without any phenomenological adjustments, which 
we refer to as realistic shell-model calculations tout court. 

The two main ingredients of realistic shell-model calculations are the A^A^ potential 
Vnn and the many-body methods for deriving Ves- As regards the latter, a main 
difficulty with its derivation is the strong short-range repulsion contained in modern 
A^A^ potentials. As is well known, the traditional way to circumvent this problem 
is the Brueckner G-matrix method. However, a few years ago a new approach to 
the renormalization of V^n has been proposed [7, 8J, which has proved to be [9] an 
advantageous alternative to the use of the G-matrix . This consists in constructing a 
smooth low-momentum potential, Viow-k, that can be used directly to derive V^s within 
the framework of a perturbative approach. 

Since its introduction in 2002, we have used this approach to derive realistic low- 
momentum effective interactions for shell-model calculations in various mass regions, 
which have provided an accurate description of nuclear structure properties. It should be 
noted that, having taken the single- nucleon energies from experiment, these calculations 
are essentially parameter- free. 

While the realistic shell-model calculations performed thus far have all been based 
on a Hamiltonian containing only two-body interactions, the role of three-nucleon 
interactions has been deeply investigated within the framework of ab initio approaches. 
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such as the GFMC and the NCSM. It is certainly a major open problem to investigate 
the contribution of three-body forces to the shell-model effective interaction. 

In this paper, we shall try to point out some questions that in our opinion remain 
open in shell-model calculations with realistic effective interactions. This is done in the 
context of a brief description of the current status of the field, which serves the purpose 
to provide the appropriate frame for our discussion. 

2. Theory and applications 

The Schrodinger equation for a nucleus with A nucleons is written as 

H^, = {H, + H,)^, = E,^„ (1) 

where 

Ho = T + U (2) 

and 

H, = Vnn - U, (3) 

T being the kinetic energy and U an auxiliary potential introduced to define a convenient 
single-particle basis. The effective interaction Ves acting only within a reduced model 
space is then defined through the eigenvalue problem 

PH,sP\^^) = P{Ho + Kfr)P|^a) = KP^a, (4) 

where the Ea and the corresponding \1/q, are a subset of the eigenvalues and 
eigenfunctions of the original Hamiltonian defined in the complete Hilbert space. The 
P operator projects onto the chosen model space, which is defined in terms of the 
eigenvectors of the unperturbed Hamiltonian Hq. 

A well-established approach to the derivation of realistic effective interactions is 
provided by the Q-hox folded-diagram expansion [9j. It is, however, outside the scope 
of this article to give any detailed discussion of this approach. Rather, our aim is to 
discuss synthetically its merits and limits. To this end, this section is divided into two 
parts. In the first one we examine in a critical way the main steps involved in the 
derivation of Vcs- In the second one we give a brief survey of some applications to show 
the accuracy of realistic shell-model calculations in the description of nuclear structure 
properties. 

2.1. Outline of theoretical framework 

2.1.1. Choice of Vnn Starting in the mid 1990s, several potential models have been 
constructed which predict almost identical deuteron observables and are phase-shift 
equivalent, namely they fit equally well (x^/datum ^ 1) the A^A^ scattering data up to 
the inelastic threshold. These are the potentials constructed by the Nijmegen group. 
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Nijm I and Nijm II [TU], the Argonne Vis potential p^, and the CD-Bonn potential 

More recently, new A^A^ potentials have been derived within the framework of 
the chiral perturbation theory following the basic idea of Weinberg [13]. Efforts in 
this direction have resulted in the construction of next-to-next-to-next-to-leading order 
(N^LO, fourth order) potentials p3l [T5] . An interesting feature of this approach is that 
three-nucleon forces arise naturally at third order (N^LO). 

We will not go further into this subject and refer the reader to for a survey of 
the current status of the chiral potentials. However, a general comment may be in order 
here. Currently there are no real efforts to go beyond the one-boson-exchange (OBE) 
model in the derivation of V^n within the traditional meson theory, the main interest 
in the field being focused on chiral potentials. As regards the latter, there remain some 
open problems, the most important ones being the three-nucleon forces beyond N^LO 
[T7] and the consistency of the Weinberg power counting [T8] . 

Of course, the existence of several different phase-shift equivalent A^A^ potentials 
raises the question of how much nuclear structure results may depend on the choice of the 
potential one starts with. However, the renormalization of the free A^A^ potential, which 
is the subject of section 2.1.2, removes to a large extent the ambiguity in the choice of 
Vat AT. We have indeed verified [9] that low- momentum shell-model effective interactions 
derived from phase-shift equivalent A^A^ potentials through the Viow-k approach do not 
lead to significantly different results. 

2.1.2. Renormalization of V^n In recent years, a renormalization-group-based 
approach has been introduced [3 E] allowing one to handle the nonperturbative 
behaviour that characterizes all realistic A^A^ potential models. In this approach, the 
high-momentum modes of Vnn are integrated out down to a certain cutoff momentum A, 
which leads to a low-momentum potential that decouples high- and low-energy degrees 
of freedom while preserving the on-shell properties of the original Vnn- By varying the 
cutoff, an entire class of low-momentum potentials can be associated with the initial 
Hamiltonian, which are all by construction softer than Vattv- This greatly simplifies, 
or even makes feasible in some cases, nuclear structure calculations. It is a remarkable 
feature of this approach that with a cutoff in the range A ~ 2.1 fm~^ low- momentum 
interactions derived from different A^A^ potential become very similar to each other. 
This has recently led to a widespread use of these low-momentum potentials in various 
contexts, going from few-body systems to nuclear matter. 

Different techniques have been developed to derive low-momentum potentials, the 
decimation of the initial potential being performed on the two-body systems in all 
cases. In our calculations, we have adopted a conventional effective interaction technique 
based on the Lee-Suzuki similarity transformation [19]. The cutoff A specifies the low- 
momentum space V within which the ^ow-k is confined and its complement Q = 1 — V. 
Then, using a similarity transformation X a new Hamiltonian = X~^HX is defined 
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which satisfies the decouphng equation 

Qnv = 0. (5) 
The low-momentum potential is then given by 

Mow-k = VHV - VTV. (6) 

The solution procedure for Eq. (5) can be found in [9]. Here we only mention that 
we adopt the iterative technique for non-degenerate model spaces proposed in [20] 
and employ a momentum-space discretization procedure with an adequate number of 
Gaussian mesh points. 

As mentioned above, the Vio„-k is developed for the two-body system, which means 
that when used in calculations for A > 2 the low-energy observables are not the same 
as those produced by the original A^A^ potential and depend, to a certain extent, on 
the value of the cutoff [21]. This dependence may be removed by complementing, for 
each value of A, the two-body l^ow-k potential with three- and higher-body components, 
which are in fact generated by the renormalization procedure when applied to an A- 
body system, even if the initial Hamiltonian contains only a two-body potential [22] . 
However, the three-body forces included in Viow-k calculations to date are generally 
approximated by fitting the parameters entering the leading chiral three-nucleon forces. 
A consistent evolution of the three-body forces in the three-nucleon system has been 
only recently achieved in the work of [23]. There, calculations are performed for A < 4 
nuclei with an initial Hamiltonian including as two- and three-body forces the N^LO 
[H] and N^LO [21] interactions, respectively. 

Note that the inclusion of forces beyond the two-body ones has been essentially 
limited to few-nucleon systems. As regards medium- and heavy-mass nuclei, only the 
A^A^ potential is used at present in the derivation of the shell-model effective interaction. 
We shall come back to this point in section 3. 

It is worth mentioning that in almost all our calculations the cutoff parameter A has 
been chosen in the vicinity of 2 fm~^. This means that l^ow~k preserves the phase shifts 
up to -^Lab ~ 350 MeV, which is the inelastic threshold. We have found [25], however, 
that the shell-model results do not change significantly for moderate variations of the 
cutoff around this value. We have also verified [25] that allowing for limited changes 
(~ 0.3 fm~^) in the value of A the ^ow-k's extracted from different phase-shift equivalent 
potentials give essentially the same results. Further study is certainly needed to better 
clarify the effect of cutoff variations on shell-model effective interactions. 

In concluding this section, it is worth insisting on the fact that the use of V|ow-k 
rather than the G matrix represents a substantial progress in the derivation of shell- 
model effective interactions. The main advantages of Vio„_k may be summarized as 
follows: i) since it is a smooth A^A^ potential, Vio„_k can be used directly in nuclear 
structure calculations within a perturbative approach; ii) it does not depend either on 
the starting energy or on the model space, as instead the case of the G matrix, which 
is defined in the nuclear medium; iii) the Viow-k's extracted from various phase-shift 
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equivalent potentials are very similar to each other, thus suggesting the realization of a 
nearly unique low-momentum potential. 

From the above, it is clear that Vio„-k is also a most valuable tool to tackle the tough 
problem of taking into account three-body forces in realistic shell-model calculations. 

2.1.3. Derivation of the effective interaction As mentioned above, the calculation of 
the effective interaction within a chosen model space P is performed by using the Q-box 
folded-diagram expansion developed for two valence particles. 
The effective interaction of Eq. (jl]) can be written as 

oo 

V,a = Q + Y.Fi, (7) 

i=l 

where the Q-hox is defined as the sum of all irreducible and valence linked diagrams 
with Viow-k replacing Vnn va. the Hi vertices, and the Fj's represent Q-box i- folded 
diagrams. The latter can be expressed in terms of Q-hox derivatives with respect to the 
energy variable a;, namely 

Fi =QiQ, 

F2 =Q2QQ + QiQiQ, 

F3 =QzQQQ + Q2QiQQ + Q2QQiQ + QiQ2QQ + QiQiQiQ, (8) 



with 

Qr. 



1 d'^Q^uj) 
ml du^ 



(9) 



ojq being the unperturbed energy of two valence particles in a degenerate model space, 
PHoP = uo. 

As a first step, one has to calculate the Q-hox. This calculation requires some 
approximations, since only diagrams up to a finite order in Hi can be included and a 
truncation of the intermediate-state summation has to be made for each Q-box diagram. 
Then, for a given Q-hox the folded-diagram series ([7]) can be summed up by the Lee- 
Suzuki iteration method [19]. 

One of the most important questions arising in this kind of approach is the 
convergence of the effective interaction expansion, to which much attention has been 
paid over the years. Before folded diagrams were taken into account, both the order- by- 
order [26] and the intermediate state [27] convergence were questioned. Furthermore, 
it was shown [28] that the presence of an intruder state in the low-lying spectrum of 
the nucleus with two valence nucleons is an essential source of divergence. The problem 
of the intermediate state convergence was substantially solved by using A^A^ potentials 
with a weak tensor force, while the introduction of folded diagrams and their summation 
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through the Lee-Suzuki scheme was proved to be an important step to handle the non- 
convergence of the order-by-order expansion of Ves, also when intruder states are present 
in the P-space [191 [29]. As a matter of fact, if we know the Q-hox and its derivatives, Ves 
is obtained after a small number of iterations. The convergence of the Q-hox, however, 
still remains to some extent an open problem, although it has been shown that the 
differences in Ke's calculated with Q-boxes to various orders, are attenuated when the 
folded diagrams are summed to all orders [30]. In particular, in [30] a comparative study 
between results obtained using a second- and third-order Q-hox in the folded- diagram 
expansion has been performed. There it was found that both calculations give very 
similar results for nuclei with two valence nucleons, which, as it was also suggested 
in [29], should be traced to the cancellation of the higher-order Q-hox diagrams by the 
folded diagrams. However, the final conclusion of [30] was that third-order calculations 
should be employed since third-order terms of the Q-box are not negligible. On the 
other hand, in almost all our calculations we have used effective interactions derived 
by including diagrams up to second order. It is worth mentioning that our choice is 
also supported by the results of [31]. There, a derived from the CD-Bonn A^A^ 

potential has been employed and the scZ-shell effective interaction has been calculated 
by summing core-polarization diagrams to all orders. The results of this calculation 
have turned out to be remarkably similar to those of a second-order one. 

There are a few other points that we would like to discuss here. The effective 
interaction constructed within the above approach contains both one- and two-body 
terms. All the one-body contributions, once summed to the eigenvalues of Hq, give 
the single-particle effective energies, which should represent the energies of the one- 
valence-nucleon system. In most of the realistic shell-model calculations to date, these 
single-particle energies are replaced by the experimental values, and only the two-body 
term of Kg is retained by employing the subtraction procedure described in [29] . 

The effective interaction is constructed for a two-valence-particle nucleus but is 
generally used also for more complex systems. This means that the three- or higher-body 
forces, which appear in the effective interaction for the latter, even if the initial potential 
is only a two-body one, are not taken into account. In p2] the lack of three-body forces, 
effective or original, was considered responsible for the failure to explain simultaneously 
binding energies and excitation spectra of several nuclei. In [33], using the two-level 
Lipkin model, it was shown that effective three-body forces are not negligible. On the 
basis of our calculations, however, we do not have at the moment a definite conclusion 
on the role of these forces. In the next section we shall comment on this point. 

Finally, it is worth recalling that the perturbation in our starting Hamiltonian is 
Hi = Viow-k — U, and, as a consequence, the Q-hox diagrams contain —U as well as 
^low-k vertices. However, in almost all calculations performed so far diagrams with the 
U potential and bubble insertions are not included. This assumption is based on the 
fact that these diagrams would cancel each other when using for U the Hartree-Fock 
self-consistent potential. The most common choice, however, is the harmonic oscillator 
potential, and the plausibility of the above approximation remains to be ascertained. 
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In some of our recent papers [HU ES] all the V|ow-k and —U diagrams up to third order 
are taken into account. 

From this discussion it seems that some questions inherent in the derivation of the 
effective interaction still remain to be settled, especially as regards the real advantages 
of an a priori more complete and consistent calculation of the effective interaction. We 
feel that further investigations aimed at finding specific answers to these questions are 
needed. However, these considerations cannot ignore the remarkable success achieved 
thus far by realistic shell-model calculations in explaining nuclear structure properties. 
This will be illustrated in section 2.2. 

2.2. Practical value of realistic shell-model calculations 

While the first efforts to derive the shell-model effective interaction from the free NN 
potential go back to more than forty years ago, the practical value of this approach 
has emerged only during the last decade or so. We shall here try to give a concise 
yet quantitative evaluation of the results produced by modern realistic shell-model 
calculations. These results concern spectroscopic properties of nuclei in various mass 
regions both near and far from stability, and have been obtained using second-order 
folded-diagram effective interactions derived from the CD-Bonn potential renormalized 
by constructing a V^ow-kwith A = 2.1 — 2.2 fm~^. In particular, calculations have 
been performed in the regions around doubly magic ^''°Sn, ^^^Sn and ^osp]^ involving 
some twenty nuclei altogether (p] and references therein, [361 137]). In practically all 
cases a very good agreement between the calculated and experimental spectra has been 
obtained. In general, the rms deviation does not exceed 100 keV. Also, the currently 
available experimental data on electromagnetic transition rates are well reproduced by 
the theory (see, for instance, [38]). 

It is well outside the scope of this article to give even a brief survey of the results 
of the realistic shell-model calculations referred to above. However, to quantitatively 
illustrate their practical value, we feel it appropriate to just discuss in some detail the 
results for the two odd-odd nuclei ^^^Sb and ^^^I. In both cases the main motivation 
for our calculations [39l [37] was to interpret new data of primary interest for the 
understanding of the shell structure of nuclei in the vicinity of ^'^^Sn. Note that ^^^Sb 
with an N/Z ratio of 1.67 is at present the most exotic open-shell nucleus beyond ^^^Sn 
for which information exists on excited states. 

In the study [39] two new transitions, in addition to the one previously known, were 
observed in ^^^Sb using 7-ray and conversion-electron spectroscopy, and a level scheme 
was constructed by comparison with the results of a realistic shell-model calculation 
performed within the theoretical framework described in the previous sections. As 
is seen in figure 1, the predicted energies are in very good agreement with measured 
values. The theory predicted also a B{E2; 6^ — )■ 4~) value which comes quite close to 
the experimental one. 

As regards the nucleus ^^''I, five high-spin excited states up to an energy of about 
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Figure 1. Experimental and calculated level scheme of ^-^^Sb. 



3 MeV were very recently identified 110] through measurements of prompt 7 rays from 
the spontaneous fission of ^^^Cf. Angular correlations were performed, but were not 
sufficient to assign spins and parities. We have interpreted the observed states on the 
basis of a realistic shell-model calculation and made spin-parity assignments accordingly 
[37] . As is seen in figure 2, a very good agreement was found between the experimental 
and calculated levels. 
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Figure 2. High-spin states in identified in [lOj compared with the calculated level 
scheme for negative-parity yrast states. 
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3. More on open problems 

In the previous sections we have already touched upon some interesting questions 
regarding realistic shell-model calculations. In this section, we give a more detailed 
discussion of what we feel are the issues worthy of further study in the near future. 

The two examples given in section 2.2 provide good evidence of the predictive 
power of realistic shell-model calculations. However, as is the case for the majority 
of our calculations, they concern nuclei neighboring a double shell closure. It is, of 
course, a very interesting problem to find out what happens with increasing number of 
valence nucleons. From this viewpoint, of special interest are long isotopic or isotonic 
chains which provide the opportunity to investigate the effect of adding neutrons or 
protons to a doubly magic core over a large number of nuclei. To this aim, we have 
conducted a study of the = 82 isotonic chain which has long been considered a 
benchmark for shell-model calculations employing effective interactions derived from free 
A^A^ potentials [HI 112] . The results of this study [9j have shown remarkable agreement 
with experiment for the low-energy spectra all over the isotonic chain. However, the 
ground-state energy, which is well reproduced for the two-valence-proton nucleus ^^^Te, 
is increasingly overestimated with increasing number of protons. It should mentioned 
that a similar behavior was found in |13] for the Sn isotopic chain. This study, however, 
differs from that of p] in two respects: the short-range repulsion of Vnn is taken care 
of by introducing the Brueckner G matrix and diagrams through third order in G are 
included in the Q-box. In [13] it was shown that the calculated binding energies turned 
out to be in good agreement with experiment by adding a small repulsive monopole 
contribution to the effective two-body interaction. 

These findings lend support to the suggestion that the overbinding of nuclei with 
more than two valence particles produced by realistic shell-model calculations may be 
traced to the lack of many-body forces, in particular the three-body ones. However, 
we have recently performed calculations for oxygen [H] and calcium [31] isotopes, as 
well as for N = 82 isotones [35], including diagrams up to third order in Vio„-k- The 
results obtained seem to indicate that third-order diagrams may play a significant role 
in curing the above "overbinding problem" . 

In summary, we may say that a deeper comprehension of the role of both genuine 
and effective three-body forces and of contributions to the Q-box beyond second order 
certainly deserves further study. In this connection, it should be mentioned that an 
attempt to explain the overbinding in neutron-rich O isotopes in terms of genuine three- 
body forces is made in [45]. 




In the calculations for ^^^Sb and ^^^I mentioned above only particle-particle and 
particle-hole matrix elements come into play. We have, however, also performed several 
calculations for nuclei with two or more valence holes [Ml H?]. While the overall 
agreement between theory and experiment may be considered satisfactory in most cases, 
it is not in general of the same quality as that obtained for nuclei in which there is no 
hole- hole interaction. At present, we don't have a clear understanding of the origin of 
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this dissimetry. 
4. Final remarks 

In this paper, we have tried to give an outhne of the present status of shell-model 
calculations employing realistic two-body effective interactions and point out some 
problems that deserve further study. To place our discussion in a proper perspective, we 
have briefly surveyed the main steps involved in this kind of calculations, whose practical 
value has been exemplified by the results of two recent calculations in the neutron-rich 
^^^Sn region. 

As regards open problems, we have focused on those issues that are raised by our 
own studies. This may have left out other interesting questions. This is certainly the 
case of the coupling with the scattering continuum which should play an important 
role in nuclei near the drip lines. We refer to |18] for a very recent study of this 
problem including references to earlier works. We hope, however, to have succeeded 
in highlighting the value of realistic shell-model calculations and the need for further 
efforts to clarify some remaining open questions. 
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